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Abstract 
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^ I The QCD sum rule for the correlator of topological charge densities xiQ'^) and 

'^ ' related to it longitudinal part of the correlator of singlet axial currents is considered 

p^ . in the framework of instanton model. The coupling constant /^/ of r/'-meson with 

^3 I the singlet axial current is determined. Its value appears to be in a good coincidence 

^^ ■ with the value determined recently from the connection of the part of proton spin 

Qs . S, carried by u, d, s quarks, with the derivative of QCD topological susceptibility 
x'(0). From the same sum rule rj — rj' mixing angle 9s is found in the framework of 

o ' two mixing angles model. The value of 9s is close to that found in the chiral effective 

ji. . theory. The correlator of topological charge densities x(Q^) at large Q^ is calculated 

D I and it was found that its Q'^-dependence matches well with Q^-dependence at low 

i-^ ' Q^ determined by the known x'(0) and by contributions of tt- and 77- mesons. 
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1 Introduction. 

Recently the vacuum expectation value of singlet axial current /q induced by external 
singlet axial field A^ has been found [1]: 

(0|j^5|0)a = 3/o%. (1) 

In ([^) j^5 is the singlet quark current 

j^x) = J2 Qixh/^^5q{x), q = u,d, s. (2) 

The term 

AL = j^,A^ (3) 



was added to QCD Lagrangian, where A^ is a constant in space and time singlet axial 
field and the limit of weak A^ field was considered. It was found in the limit of massless 
u,d, s quarks [jl], |[: 

f^ = (2.8 ± 0.7) ■ 10^2 GeV^. (4) 

This result was obtained, constructing the QCD sum rule in the external field A^ to 
determine the part of proton spin S carried by u, d, s quarks and related to the proton 
matrix element of the current j^s : 

2ms^S = {p,s\ j^5 \p,s), (5) 

where s^ is the proton spin 4- vector, m is the proton mass. The sum rule for S essentially 
depends on /q and the numerical result (4) comes in two ways, resulting in the same 
value: 

1) from the requirement of consistency of phenomenological and calculated in QCD 
sides of the sum rule as functions of the Borel parameter M^; 

2) by using for E its experimental value, S = 0.3 ± 0.1. 

As was shown in |l[| , in the limit of massless u, d, s quarks /q is related to the first 
derivative x'(0) of the correlator of topological charge densities Q^i^x) 

X{q') = I Jd'xe^'^^O I T{Q,{x), Q,{0)} \ 0) , (6) 



a 



Q^i^) = ^G'" (x)G" (x), (7) 



in 



where G^^,{x) is gluonic field strength, GJ]^(x) is its dual, G"^^ = {l/2)e^y\^G'^^\ 

fo = 12x'(0) . (8) 

As follows from (^), 

X'(0) = (2.3 ± 0.6) • 10-^ Gel/^ _ ^g) 



Let us remind the derivation of (^. Using (j^) we can write: 

(0 I j^5 I 0)^ = lim t fd^xe^'^O \ T{j,,{x), j^^iO)} I 0)A, = \imP^,{q)A,. (10) 

q^O J q-*0 

The general structure of P^v{q) is 

P^M = -PL{q')S,. + PT{q'){-S,uq' + q^Qu). (11) 

Because of anomaly there are no massless states in the spectrum of singlet polarization 
operator P^^, even for massless quarks. Pt',l(q'^) also have no kinematical singularities at 
g^ = 0. Therefore, the nonvanishing value Pfj,u{^) comes entirely from PL{q^). Multiplying 
P^iv{.q) by q^qu and using the anomaly condition 

d^j^^{x) = 2NjQ^{x) + 2iY^ m,q{x)^^q{x) (12) 

{Nf is the number of flavours, Nf = 3), in the limit of massless quarks we get: 

q,quP,u{(i) = -PL{q')q' = 36x(g'). (13) 

As is known 0, x(0) = if there is at least one massless quark. (H) follows directly from 
(P, (|10D, (In]), (|i|D. According to (|ID we have also: 

/o = -(1/3)Pl(0). (14) 

The attempt to determine /q directly by constructing a special QCD sum rule for this 
aim was performed in |^. However, this attempt failed: it was found that the operator 
product expansion (OPE) used in the sum rule does not converge on the accounted in 
calculation terms - the unaccounted higher order terms of OPE should be of importance. 
In the present paper we use the instanton model (as review see 0) for calculation 
of these higher order terms. The idea that instantons give the main contribution to the 
longitudinal part PL{q^) of the correlator of singlet axial current P^uiq) at intermediate 
I g^ |~ 1 GeV^ is not new - it was suggested as early as in 1979 [0, §. In it was 
argued that the appearance of r]- and T^'-mesons as almost pure octet and singlet states in 
SU{3) flavour symmetry, i.e. the large mixing of uu + dd and ss in this channel, cannot 
be described by perturbative QCD and may be attributed only to dominating instanton 
contribution. 

The sum rule for direct determination of /q is constructed. At standard parameters 
of the instanton model the value of /g found from the sum rule is in a good agreement 
with (I). 

For massless quarks the phenomenological side of the sum rule is saturated by r/'- 
meson contribution (plus contributions of excited states, approximated by continuum). 
The strange quark mass rus may be also accounted in the sum rule. In this case the 
contribution of 77-meson also comes into play and rj — t]' mixing angle can be found from 
the analysis. In the model of two mixing angles the value of the largest angle 9s is 
determined and appears to be in a good agreement with the values found from the chiral 
theory and phenomenology 
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In the framework of the same instanton model the g^-dependence of the topological 
charge densities correlator xiQ"^) is determined at space-like Q^ = —q^ > 0. At inter- 
mediate Q^ ~ IGeV"^ it matches well with the curve of the xiQ"^) behaviour at low Q^, 
found in |]12[ on the basis of (^ and contributions of Goldstone bosons vr" and 77. 

2 The sum rule. 

Strictly speaking, the quantity /g given by (|l]) is defined as nonperturbative part of the 
induced by the external field A^ vacuum expectation value (0 | j^s | 0)a with perturbative 
contribution subtracted. The numerical value (^ corresponds just to such definition. The 
reason for this definition is that in the sum rule for S, from which the value of /g (|^) 
was determined, all perturbative contributions were accounted explicitly and only non- 
perturbative part of (0 I j^5 I 0)^ was parameterized by unknown constant /g. Similarly, 
x(g^) in (P) and (|^) has the meaning of nonperturbative part of the correlator of topo- 
logical charge densities. Such definition is physically reasonable, since the perturbative 
part of @ is badly divergent, strongly depends on renormalization scheme, and therefore 
has no physical meaning. The same statement refers to perturbative contribution to /g. 
Such separation of perturbative and nonperturbative contributions allows one to avoid 
any uncertainties in the sum rule for the physically measurable quantity S. 

The idea of determination of /g or, what is equivalent, proportional to /g quantity 
-Pl(O) was suggested in 0]. In short, it was the following. 

The imaginary part of Piio^) is represented by contributions of the lowest resonance 
- r^'-meson - and continuum: 

ImPL{q^) = 3nf^,ml6{q^ - mj,) + Piq^iq' - So) . (15) 

Here /^, is the coupling of r^'-meson with the singlet axial current 

(0 I jf,5 I v') = ^v^/v^M (16) 

in the limit of massless u, d, s quarks (g^ is 77' momentum). The second term in the right 
hand side (rhs) of (15) represents the contribution of continuum, sq is the continuum 
threshold. The continuum contribution corresponds to gluonic bare loop in the correlator 
(^ and is equal to 

/3tf) = g.^ (17) 

In order to get the nonperturbative part of IuiPl the perturbative part equal to /?(g^) 
must be subtracted from (15) what gives: 

ImPLiqXonp = 37rf^,ml,S{q^ - mj,) - P{q^)e{so - g') . (18) 

As was shown in [^, the nonperturbative part of -Pl(O) is given by 

Pmnou, = -SP^> + ll ^ds . (19) 





Therefore the problem reduces to determination of the couphng constant /3. This can 
be done by standard technique of the QCD sum rule approach. Let us write OPE for 
Pl{Q^) = (36/(5^)x(Q^) at large Q^. We account the instanton contribution in the 
instanton liquid approximation. (The instanton contribution was not accounted in 0.) 
The OPE for PLiQ^) has the form [g: 
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+ ^73—.9r"'{0 I Gl^C^.G',^ I 0)+ 
9a? 



I "'-^s fabc rade In \ /~fb r-ic r^d r^e i -\c\r^h r^c r^d r-ie \ r\\ i 
2~06 ^ ' ^'^ a/3 A^i^ "/3 ^ ^^^ iia^ av'^ ^ifi^ iiv I U; + 

+ 18g2 I dpn{p)p'Kl{Qp) . (20) 

In (^) operators up to dimension 8 were accounted. The first term in the rhs of (^) is 
the bare loop contribution, the last one is the contribution of instanton |]^, Q, K2{x) is 
the McDonald function, p is the instanton size and n{p) is the instanton density. (In our 
normalization anti-instantons are accounted in the coefficient 18 in front of the last term 



in (pOD and n{p) represents the instanton density.) For n{p) we use the Shuryak model 



[^, ^] of instanton density: 

n{p) = no5{p - pc) . (21) 

As was demonstrated by Shuryak and his collaborators 0, this model well describes 
many hadronic correlators in QCD. For numerical values of parameters in ( pi]) we choose: 
riQ = 0.75 ■ 10~^ GeV^, pc = 1.5 GeV^^, inside allowed by this model limits. Particularly, 
at this no the standard value of gluonic condensate 

rv 
(0 I —G^ I 0) = 0.012 GeV^ (22) 

vr 

may be attributed entirely to instantons. In order to estimate 8-dimensional gluonic con- 
densate we assume the factorization hypothesis - the saturation by vacuum intermediate 
states. Then p: 

/abc£ade/r\ \ /^b r^c /^d r^e i i n/^^ r^^ r^^ r^^ I C\\ 
/ \U I '^fj.u'^al3^fiu^al3 + ^^^ fia^ au^ fj./d'^ /3u I ^/ " 

= ^{0 \ g;,g;^ \ 0)' . (23) 

It should be mentioned that the calculation of the same term in the instanton model 
would give quite different result: 

(2^1 • 37^/7)r^op,-^ (24) 

which is by an order of magnitude larger than ( PBj ) at accepted model parameters. This 
fact is not surprising. Indeed, for the gluonic condensate with k gluonic fields on dimen- 
sional ground we would have: 



{0\G'\0)r^ Jn{p)^j-^dp, (25) 

and the integral (^) diverges at high enough k at any physical n{p). Therefore, one may 
expect that the instanton model overestimates the value of 8-dimensional gluonic conden- 
date and accept the estimation (p3D based on factorization hypothesis. This estimation 
is supported by the analysis of the sum rules for heavy quarkonia with the account of 
8-dimensional operators |T3| where the factorization hypothesis was used. Much larger 
values of 8-dimensional gluonic condensates would contradict the analysis in |T^. Of 



course, two or, may be, even 3 times larger values as ( P5| ) are not excluded, but, luckily, 
the contribution of these condensates to the sum rule is small and even its 3 times in- 
creasing does not influence seriously the result. For the 6-dimensional gluonic condensate 
there is no other independent estimation as given by instanton model [§|: 

^' ^r'^iO I Gl^GlC^^^ I 0) = ^1(0 I ^G' I 0) . (26) 



The phenomenological representation of Pl{Q^) follows from (p!5|). Equating these two 
representation and applying the Borel transformation to both sides of this equality we get 
the sum rule: 

p2 2 -vn?i/M'^ 




+ ^^(0 I ^G^ I 0)^ + \%n,p\Bu^Q^Kl{Qp:), (27) 

Ei(x) = l-(l + a;)e-", (28) 

where Bm'^ means Borel transform. In ( P?] ) e corresponds to contribution of 6-dimensional 
gluonic condensate. The instanton model estimation (|26|) gives e = 2.2 GeV^. Since we 
expect that the instanton model overestimates 6-dimensional gluonic condensate also, we 
put e = 1 GeV^ and include possible uncertainty in the error. The Borel transformation 
of McDonald function can be done by using its asymptotic expansion, what leads to (see 



H): 



18noptB,r.Q'Kl{Qp,) ^ 9n,M'plV^e~'^'^'^ (m'pI + ^ + ^^j • (29) 

In our M^ domain the next terms of the expansion are small. In order to verify this fact 
numerical calculation was done, using integral representation of McDonald function. 

The results of the calculation of f^, according to the sum rule ( pTf ) are plotted in Fig.l. 
(It was put Aqcd = 200 MeV, a^ was calculated in the leading order approximation.) 
The standard estimation procedure of the M^ interval where the sum rule is reliable - the 
requirement that highest order terms of OPE are small - does not work here, because the 
instanton contribution dominates, it comprises about 75-80% of the total. So we go to 
physical arguments. We put the continuum threshold sq = 2.5 GeV^ close to the position 
of the second resonance with rj' quantum numbers, rj' {lAAO) (probably, r7(1295) the belongs 
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to octet), and require that the continuum contribution to bare loop does not exceed 
~ 50%. As a low limit of M^ interval we choose the M^ value where the M^-dependence 
starts to rise steeply. These requirements result in 1.2 < M^ ^1.6 GeV^. In this interval 
M^-dependence is not very strong and we estimate /^, ~ (2.4 ± 0.6) ■ 10"^ GeV^. (The 
error includes 15% possible variation of pc-) The contribution of the second term in the 
rhs of (|19D is negligibly small. From ([1^) and (^ we have finally: 

f^ = (2.4 ± 0.6) ■ 10^2 GeV^ , (30) 

X'(0) = (2.0 ± 0.5) • 10"^ Ge1/^ (31) 

in a good agreement with @ and x'(0) value found in |[l|], x'(0) = (2.3 ±0.6) • 10"^ GeV"^. 

3 The account of strange quark mass, rj' — rj mixing 
angle. 

Consider the polarization operator Piiff) with the account of the strange quark mass mg 
and determine the coupling constant /^/ of physical tj' . The u- and d-quark masses are 
disregarded as before. Using the definition of P^v{q) and (|12D we have: 

-PL{q')q' = q,q.P,M = i J d^xe^'^^ {0 \ T{2NfQ,{x), 2NfQ,{0)}+ 
+T{2NfQ,{x), 2im,s(0)75s(0)} + T{2zm,s(x)75s(x), 2NfQ5{0)}- 

- 4m^T{s(x)75s(x), 5(0)755(0)} | 0) + 4m,(0 | s(0)s(0) | 0) . (32) 



The last term in ( |5^ ) is caused by the equal-time commutator. Perform the OPE in the 
rhs of (0). Then in comparison with (|20|) three additional terms appear: the equal-time 
commutator term, proportional to m^ term, corresponding to bare loop of strange quarks, 
and arising from the second and third terms in ( ^2]) the term, proportional to quark-gluon 
condensate [|l5ll: 



^(0 I sa,,iX^/2)G;,s I 0) = ml{0 \ ss \ 0), (33) 

GeV^ was determined in 
sum rule is now: 



where ttt-q = 0.8 GeV^ was determined in |jT^. After Borel transformation the rhs of the 



RiM^)-AmAO I .-. I 0) + §M^^o(^)i.-«/^ -6a.^(0 | ss | 0), (34) 

where R(M^) us the rhs of (^7|) and 

Eo{x) = 1 - e-^ (35) 



L = /n(MVA')//ri(/iVA') . (36) 

The factor L~®/^ accounts quark mass anomalous dimension. 

It is useful to consider also the correlator P^uil) in the case when one of the currents 
is still j^5{x), but the other is the current of u- and tZ-quarks: U'^fTl^u + d'j^'-f^d. In this 
case the rhs of the sum rule is equal to 

lR{M^)-2a.^{0\rss\0). (37) 

In the phenomenological - left hand side (Ihs) of the sum rule - rj'- and r^-mesons are 
contributing and their mixing must be accounted. We adopt the two mixing angles model 
[^ , which is based on the low energy chiral effective theory and describes the experimental 
data better 0, |lOl than the one mixing angle model. In this model the couplings of rj- 
and 77'-mesons /,, and /,,' to octet and singlet axial currents are related to the couplings 
of fictitious octet and singlet pseudoscalar states /§ and /i by: 

f^ = fsCoses f^ = -f,szn9, 

f^, = fssines f^ = f,cosd, (38) 

The 77'- and ry-meson contributions to PiiQ"^) can be easily calculated in this model. It is 
convenient to present them separately for the cases when one of the currents is 57^755 or 
U'^y.l^u + dj^'j^d (the other is always j^s). Instead of the Ihs of (p7|) we have now: 
for S'j^'j^s current: 



2 -m^,/M^ 



m:,e """"Vif 005^6*1 - V2{fs/ fi)sin98Cos9i ] + 



+ mle-</^'' f^[ szn% + V2{fs/f\)stne,coses ]; (39) 

for ^7^75^ + d'jfj^'j^d current: 

Vi [ 2cos% + V2{fs/fi)stnescosei ] + 



2 -m?,/NP 

m„,e ^ 



+ mle-<'^''"fl [ IsinH^ - V2{fs/ fi)cos9s smd, ] . (40) 

Taking the sum of (p9D, (^0]), putting 9i = 9^ = and equating it to (|3^) at m^ = 0, we 
get the previous result with /i = /q. At nonzero rris the mixing angles must be accounted 
and for the sum of (^9]), (EO) we get: 



3m^,/i 



2 r2 -m^ilM' 



2 /J\,f2 



2 

COS 9i H 7^e- ^ ^" sin 9 1 



3mi 9 f So \ o/Q m^irti 



R{M') - 4m,(0 I ss I 0) + i^^ M-^)^'^' - ^^-^^^^(O \ss\^). (41) 



Take now the difference of (153) and one half of (HDl). The corresponding sum rule is: 
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-^m;,e-</^Vi/8 



2 



sinOs cosOi ^e -j' '"'' cosOg sinOi 



m\ 



_4,„,(0 I ,-, I 0) + gAf^E,(-^]L-/» - 3a.!^(0 | ... | 0) . (42) 



The theoretical value of Q\ found in [^, [1^, |TT[ is small: Q\ = —(2.7° — 4°). (The phe- 
nomenological value |TD|, |TT[] is slightly higher: 6i = —9.2°). Therefore, with a good 
accuracy we can put ^i ^ in (p|),(^2D. Then (^) determines /f ~ f^,. M^-dependence 
of ff is presented on Fig.l. (The numerical values (0 | ss | 0) = —1.11 ■ 10~^ GeV^ and 
ms{l GeV) = 150 MeV were used). From the curve in Fig.l the estimation follows: 

f^, = (3.2 ± 0.6) ■ 10^2 GeV\ f^> = 178 ± 17 MeV. (43) 

The ratio of (^2|) to ( PD gives the value of mixing angle Og. In the approximation 6*1 ~ 



and at fs/fi = 1.12 g |T0|, |TT1 it is equal: 

^8 = -(17.0 ±5.0)°. (44) 

The account of ^i = —2.7° changes 6*8 to 

^8 = -(18.8 ±5.0)°. (45) 

The values of /^/ and 6^ mixing angle are in an agreement with ones found in [^ |TU|, |TT| 
from the low energy effective theory or phenomenology. They are correspondingly: 

theory: /^/ = 151 MeV; 9s = -21° . 

phenom.: /^/ = 153 MeV; 0^ = -21° . (46) 

4 Q^-dependence of xiQ'^)- 

Using OPE for PiiQ"^) - (0) and (|g), the Q^.dependence of x{Q^) at high Q^ can be 
found. Since, by definition of x((5^), the perturbative part should be subtracted, the first 
term in the rhs of (^) is omitted and we have: 

xiQ'') - -^(0 I f G^ I 0>(l + ^) - ^.a.i,(0 I ^G^ I 0>^- 

- inoQV'AKOft), (47) 



where e parameterizes the 6-dimensional gluonic condensate contribution (see (pT])) and 
( p3D was used. x(Q^) (0) ^^ plotted in Fig. 2. It is instructive to compare x{Q^) at large 
g2 with xiQ^) at low g2 found jn |12|: 



xm = x(0) - x'io)Q' - -^flQ' 

where 



'niu-mdV ml 1 m^ 



+ ::- 



jriu + mdj Q'^ + ml 3Q'^ + m'^_ 



(48) 



X(0) = ^M^(o I uu I 0) (49) 

rriu + md 

and the last term represents the contributions of 7r°- and ry-mesons. The curve of the low 
Q^ behaviour of x{Q^) (0) is also plotted in Fig. 2, for x'(0) i^ ^^s chosen the value found 
in |jl|: x'(0) = 2.3 ■ 10"'^ GeV^. As is seen from Fig.2, both curves matches rather well in 
the domain Q^ ^ 0.4 - 1 GeV'^. 

5 Discussion. Comparison with results of other 
works. 



As was mentioned in the Introduction, instantons are the most plausible QCD objects 
for description of physical r]' channels, or, what is equivalent, of the longitudinal part 
of singlet axial vector current correlator. The calculations of the pseudoscalar current 
Js = ^750' correlator with rj' quark content performed by Shuryak and Verbaarschot \T7\ 



and Schafer [|I^ in the framework of various instanton models demonstrated that this 
correlator is described in an agreement with its phenomenological coordinate dependence 
up to distances x ~ 0.3 fermi (and in some of them even up to larger ones). Therefore, 
beforehand, we could expect that the instanton model is suitable for consideration of the 
problem in view. However, in this paper we used the simplest version of the instanton 
model - the instanton liquid approximation with instanton density given by the Shuryak 
model (^). For this reason the accuracy of our results is limited. 

Since the main contribution (about 80 %) to the sum rules from which /q and /^, were 
found comes from the instanton term, the main uncertainty is caused by the instanton 
parameters no and pc- These parameters were taken from the best fit to various hadronic 
correlators, as well as some other QCD objects, like gluonic condensate performed by 
Shuryak and his collaborators [§ . Possible uncertainties are included into the errors. The 
errors in the determination of the mixing angle 6s are smaller, because instantons do not 
contribute to the rhs of (42). If for /i and /s their phenomenological values will be taken: 
/i = 1.28/7r, fs = 1-15/^ [^, |TD|, |TT|] instead of /i = /^/ found from the sum rule, then we 



would have from (||) (at 9i = -2.7°): 

^8 = -(26.5 ±3.5)°. (50) 

The obtained above value of /^/ ( |i3| ) is a bit higher than the low energy data (|46D, but, 
taking in mind the uncertainties, one may consider the agreement as satisfactory. 

The a^-corrections were not accounted in our calculation. They are essential in the 



case of the first term in the rhs of ( pT]) [19 , but this term contributes only 5% to the total 
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result. The as-corrections to gluonic condensate contributions are masked by uncertainties 
in higher order gluonic condensates. Among the terms proportional to rris and ml, the 
a^-corrections appear only to the two last terms in (0), not to the main, proportional to 
rris, term — 4771^(0 | ss | 0). They do not much influence the value of the mixing angle. 
There are also instanton corrections to the term ~ m^ of the same order of magnitude. 
The account of all these corrections to the ~ m^ terms is the subject of further study. 
We believe that after their account the value of the mixing angle will be still in the limits 
of accepted above errors. 

The slopes in M^ of the left and right hand sides of ( plD are different - positive in 
the Ihs and negative in the rhs. For this reason it is impossible to determine m'i, by 
differentiation of (^T]) in M^ as it was done sometimes in the QCD sum rule approach: 
whereas the sum rule is satisfactory, its derivative is not. 

In their recent paper Narison, Shore and Veneziano |T^ studied the problem of deter- 
mination of x'(0) or /^, in the framework of the standard QCD sum rule approach with no 
account of instantons. Their result for x'(0) or /3 is essentially - 3-4 times - smaller than 
ours. The principal difference from the presented here calculations (besides instanton 
contribution) is that the authors of |TP[ have chosen much larger values of the continuum 
threshold sq and of the effective Borel parameters: sq = 6GeV^ and M^ ~ 3 — AGeV"^. 
Therefore, they accepted the model of hadronic spectrum in J^*^ = ^ flavour singlet 
channel with a gap between ?7'-meson mass and 2.5 GeV. However, there are at least 
3 resonances with t]' quantum numbers between r]' mass and 2.5 GeV. That is why we 
think that such a model is not acceptable physically. The other drawback of the sum rule 
used in |T^ (eq.(D.ll) of |jl9|, which is similar to (0), but without rj — rj' mixing) is that 



the main contribution to the sum rule comes from the terms proportional to rris and ml. 
These terms comprise 60% of the flnal answer. This means that SU(3) flavour symmetry 
is badly violated in the sum rule |T^ in contradiction with experiment. Moreover, if we 



introduce rj — rj' mixing (what was not done in ||T9[), we may calculate the rj — rj' mixing 
angle, representing the phenomenological side of the sum rule by (|5^),(H0). Then the 



result for the mixing angle Og is: 6s ~ 45 , i.e. rj' is not mainly flavour singlet and rj is 
not mainly octet - in evident contradiction with experiment. 

In conclusion, we have shown that the instanton model even in its simplest version 
describes reasonably well the properties of the topological density correlator - x'(0) and 
x(Q^) at large Q^, the values of the r]' coupling constant and the r]' — rj mixing angle. 

This work was supported in part by RFBR grant 97-02-16131. 
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Figure captions. 



Figure 1. Functions f'^i{M'^), /^,(M^) determined by equations (27), (41) correspond- 
ingly. In (41) Oi = 0. 

Figure 2. Function x(Q^) ^^ low and large Q^ - solid lines (see equations (48) and 
(47) correspondingly). The dashed line represents matching curve drawn by hand. 
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